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Abstract. Plane-symmetric gravitational waves are considered as gravitational 
lenses. Numbers of images, frequency shifts, mutual angles, and image distortion 
parameters are computed exactly in essentially all non-singular plane wave 
spacetimes. For a fixed observation event in a particular plane wave spacetime, the 
number of images is found to be the same for almost every source. This number 
can be any positive integer, including infinity. Wavepackets of finite width are 
discussed in detail as well as waves which maintain a constant amplitude for all 
time. Short wavepackets are found to generically produce up to two images of 
each source which appear (separately) only some time after the wave has passed. 
They are initially infinitely bright, infinitely blueshifted images of the infinitely 
distant past. Later, these images become dim and acquire a rapidly-increasing 
redshift. For sufficiently weak wavepackets, one such "flash" almost always exists. 
The appearance of a second flash requires that the Ricci tensor inside the wave 
exceed a certain threshold. This might occur if a gravitational plane wave is 
sourced by, e.g., a sufficiently strong electromagnetic plane wave. 



1. Introduction 

The theory of gravitational lensing has by now reached a considerable degree of 
sophistication [1, 2, 3, 4]. Theorems have been found predicting (or bounding) the 
number of images in very general systems [5, 6, 7]. Shapes of stable caustics have 
been exhaustively classified [2, 3, 8, 9], a non-perturbative notion of the lens map has 
been obtained [10], and various universal behaviors of images have been found for 
sources lying near caustics [11, 12, 13]. These general results have been complemented 
by a number of detailed calculations for specific types of lenses. The majority 
of such calculations have been performed within the quasi-Newtonian viewpoint of 
gravitational lensing commonly used in astrophysics [3]. While various fiavors of this 
formalism exist, most require that bending angles be small and that all lenses be 
nearly-Newtonian mass distributions. 

This is to be contrasted with the more fundamental picture of gravitational lensing 
where light rays are modelled as null geodesies in a Lorentzian spacetime. Within 
this context, specific lensing calculations have been performed in Kerr, Reissner- 
Nordstrom, and a handful of other geometries [1]. While curvatures and bending 
angles may be large in these examples, they all involve (at least conformally) stationary 
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spacetimes. It is of interest to understand if qualitatively new effects appear in 
dynamical cases. 

This paper considers the bending of light by (decidedly non-stationary) plane- 
symmetric gravitational waves. Gravitational lensing by gravitational waves has 
previously been considered by a number of authors, although almost all of this work has 
been carried out within the weak-field regime [14, 15, 16, 17, 18]. One exception is [19], 
where redshifts were computed in an exact solution to the vacuum Einstein equation 
representing a plane gravitational wave. This work considered only very specific 
waveforms, and was confined to a coordinate patch too small to include caustics and 
other effects associated with the formation of multiple images. Separately, extensive 
work has been devoted to non-perturbatively understanding the geodesic structure of 
generic plane wave spacetimes [20, 21, 22, 23, 24, 25, 26]. This is clearly a subject 
closely related to gravitational lensing, although few explicit relations between the two 
subjects appear to have been made (see, however, remarks in [1]). 

It is the purpose of this work to provide a comprehensive and non-perturbative 
discussion of lensing in plane wave spacetimes. These geometries are a well-known 
subclass of pp- waves; plane-fronted waves with parallel rays. Many plane wave 
spacetimes are exact solutions to the vacuum Einstein equation. Others may be 
interpreted as, e.g., exact solutions in Einstein-Maxwell theory. While rather idealized 
from the astrophysical perspective, plane wave spacetimes admit a wide variety of 
interesting phenomena. Depending on the waveform, any number of astigmatic and 
anastigmatic caustics may exist. Examples admitting any specified number of images - 
even an infinite number - are easily constructed. In particular, even numbers of images 
can exist [1]. The image count can also change in time even when a source does not 
cross an observer's caustic. Despite all of these properties, plane wave spacetimes are 
geometrically very mild. They are topologically equivalent to and admit coordinate 
systems which cover the entire manifold. 

Aside from their value as models of gravitational radiation, the plane wave 
spacetimes considered here have also found numerous applications via the Penrose 
limit. This limit provides a sense in which the metric near any null geodesic in 
any spacetime is equivalent to the metric of an appropriate plane wave spacetime 
[28, 29]. It allows problems in relatively complicated spacetimes to be reduced to 
equivalent problems in plane wave spacetimes (which are often much simpler). This 
has been particularly valuable within string theory and related fields [30, 31]. Penrose 
limits have also been applied to ordinary quantum field theory in order to investigate 
causality and effective indices of refraction for photons and gravitons propagating in 
curved spacetimes [24, 25, 32]. More recently, Penrose limits were used to deduce 
the effect of caustics on Green functions associated with the propagation of classical 
fields in curved spacetimes [23]. Given the content of the Penrose limit, lensing in 
plane wave spacetimes might be related to lensing in generic spacetimes as seen by 
ultrarelativistic observers. We make no attempt to justify this conjecture, however. 

This paper starts by providing a self-contained review of plane wave spacetimes 
in Sect. 2. Although most of this material is not new [20, 21, 22, 23, 24, 25, 26, 27], 
it is not widely known. Sect. 3 then derives the number of images of a point source 
that may be viewed in plane wave spacetimes. Under generic conditions, this is found 
to depend only on the waveform and a certain time parameter associated with the 
observation event. The number of images does not depend on any properties of the 
source. Once this is established. Sect. 4 computes image positions, frequency shifts, 
angles, and image distortion parameters in general plane wave spacetimes. Sect. 5 
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applies these results to symmetric plane waves, which have constant waveforms. These 
geometries produce an infinite number of images of almost every source. Their lensing 
properties are found to change significantly if the Ricci tensor is increased beyond a 
certain threshold. Lastly, Sect. 6 discusses "sandwich waves;" wavepackets with finite 
width. These spacetimes generically admit images which appear at discrete times and 
then persist indefinitely. Such images initially provide infinitely blueshifted, infinitely 
bright pictures into the infinitely distant past. Very quickly, however, such images 
become highly redshifted and effectively fade away. One of these "transient fiashes" 
is produced by almost every sufficiently weak vacuum (Ricci-ffat) wave. For weak 
waves, a second ffash appears only if the Ricci tensor of the wavepacket exceeds a 
certain threshold. Throughout this work, the spacetime is assumed to be everywhere 
transparent. The language used also assumes that the geometric optics approximation 
[3] holds even in situations where it would be severely strained (such as when light is 
emitted near an observer's caustic). 

Notation 

This paper restricts attention only to plane wave spacetimes in four spacetime 
dimensions. Our sign conventions follow those of Wald [33]. The signature is — h++. 
Latin letters a, 6, . . . (and occasionally ^4, 5, . . .) from the beginning of the alphabet 
are used to denote abstract indices. Greek letters /i, i^, . . . are used to denote four- 
dimensional coordinate indices. The Latin letters j, . . . are instead coordinate indices 
associated with the two directions transverse to the direction of wave propagation. 
Objects involving the latter type of index are often written in boldface with all indices 
suppressed. They are then manipulated using the standard notation of linear algebra 
[e.g., AkiBkj = (A'''B)^j and |x| = ^yxiXi]. Overall, notation related to plane wave 
spacetimes closely follows the conventions of [23]. 

2. Geometry of plane wave spacetimes 

Plane wave geometries may be interpreted as simple models for gravitational waves 
emitted from distant sources. Alternatively, they arise as universal limits for the 
geometries near null geodesies in any spacetime [28, 29]. The typical definition of 
a plane wave spacetime (M, ^^b) requires that M = and that there exist global 
coordinates {u^v^x.) = {u^v^x^ ^x'^) : M ^ such that the line element takes the 
form 

g^^dx^dx"" = -2dudv + Hij{u)x'x^du^ + {dx^f + {dx^f. (1) 

Hij = is any symmetric 2x2 matrix. Its components describe the waveforms 

associated with a wave's three polarization states^. The u coordinate is interpreted as 
a phase parameter for the wave, while v affinely parametrizes its rays. The remaining 
two coordinates span spacelike wavefronts transverse to the wave's direction of 
propagation. 

J One of the three polarization states associated with plane wave spacetimes vanishes in the vacuum 
case Rab = 0. This leaves the usual two polarizations associated with vacuum general relativity. 
Note that the six polarization states typically stated to exist for linearized "plane" gravitational 
waves in generic theories of gravity [34] cannot all be represented by the metric (1). Three of these 
polarizations may be realized only as geometries which are rather less plane-symmetric than those 
considered here. 
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Note that if H = in some region, the spacetime is locahy flat there. In terms 
of a Minkowski coordinate system (t and V satisfy u = {t ^ x^) / \/2 and 

V = (t — x'^)/\/2 in such a region. We consider only nontrivial plane waves, so H 
cannot vanish everywhere. 

The physical interpretation of the so-called Brinkmann metric (1) as a plane- 
symmetric gravitational wave follows from considering the integral curves of the 
vector field ^^ = {d/dvY. These curves form a null geodesic congruence which 
may be interpreted as the rays of the gravitational wave. Va^^ = 0, so these rays 
have vanishing expansion, shear, and twist. There is therefore a sense in which 
they are everywhere parallel to one another. All rays are also orthogonal to the 
family of spacelike 2-surfaces generated by the two commuting spacelike vector fields 
■^{i) ~ {d/dx^Y. The induced metric on each such surface is fiat: The wavefronts 
are 2-planes. The curvature is constant on these planes in the sense that the X^^-^ are 
curvature collineations: 



^X^i^Rabc^ = 0. (2) 

Despite this, the X^.^ are not everywhere Killing. There do, however, exist linear 
combinations of X^.^^ and £^ which are Killing, itself is also Killing, which may be 
interpreted as a statement that plane waves do not deform along their characteristics. 
Plane wave spacetimes admit a minimum of five linearly independent Killing fields. 
Note that in fiat spacetime, five (out of the total of ten) Killing fields are symmetries 
of all electromagnetic plane waves [27]. Killing fields of plane wave spacetimes are 
discussed more fully in Sect. 2.3. 

All non-vanishing coordinate components of the Riemann tensor may be 
determined from 

Ruiuj ^ij' (^) 

It follows from this that the Ricci tensor is 

i?a5 = -TrH44, (4) 

where Tr denotes the ordinary (Euclidean) trace of the 2x2 matrix H. The Ricci 
scalar always vanishes in plane wave spacetimes. More generally, there are no nonzero 
scalars formed by local contractions of the metric, the curvature, and its derivatives: 
R'^Rah = RabcdR''^'''^ = Rabcde^'^^^Rfh'''^ = . . . = 0. This is analogous to the fact that 
plane electromagnetic waves satisfy, e.g., FabF^^ = e^^'^^ Fab Fed — 0- Note, however, 
that electromagnetic plane waves are not the only electromagnetic fields with vanishing 
field scalars. Similarly, plane wave spacetimes are not the only curved geometries with 
vanishing curvature scalars [35]. 

It follows from (4) that plane wave spacetimes satisfying the vacuum Einstein 
equation (and the vacuum equations of many alternative theories of gravity [36]) are 
characterized by the simple algebraic constraint TrH = 0. For vacuum waves, there 
exist two scalar functions /i+ and /ix such that 

hj^ and /ix describe the waveforms for the two polarization states of a gravitational 
plane wave propagating in vacuum. A plane wave is said to be linearly polarized if /i+ 
and /ix are linearly dependent (in which case one of these functions can be eliminated 
by a suitable rotation of the transverse coordinates x). 
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If hj^ and hx have compact support, the geometry is said to be a sandwich wave. 
This name evinces the image of a curved region of spacetime "sandwiched" between 
nuh hyperplanes in a geometry that is otherwise Minkowski. Physicahy, it corresponds 
to a wavepacket of finite length. Note that the planar symmetry considered here is 
very special in the sense that passing waves do not necessarily leave any "tail" behind 
them. After interacting with a sandwich wave, all observers enter a region of spacetime 
which is perfectly fiat. There is a sense in which test fields propagating on plane wave 
spacetimes also have no tails [23, 37]. 

A general (not necessarily vacuum) wave profile H may be built by adding to (5) 
a term proportional to the identity matrix 5. There then exists a third polarization 
function h\\ such that 

H= f y 

If the Ricci tensor of such a wave is associated with a stress-energy tensor via Einstein's 
equation, that stress-energy tensor obeys the weak energy condition if and only if 
h\\ > 0. Assuming this, the stress-energy tensors associated with (6) are very simple. 
They could be generated by, e.g., electromagnetic plane waves with the form 

Fab = 2h\v^auV^x\ (7) 

Alternatively, (6) could be associated with the stress-energy tensor of the massless 
Klein- Gordon plane wave 

0= / h\{w)dw. (8) 

J u 

Besides the vacuum case h\\ = 0, another interesting class of wave profiles are 
those that are conformally flat. These satisfy /i+ = /ix = 0, so H cx 5. As 
gravitational lenses, all caustics of conformally-flat plane waves are associated with 
"perfect" anastigmatic focusing. For more general plane waves, caustics are typically 
(but not necessarily) associated with astigmatic focusing. 



2.1. The matrices A and B 

The geometry of plane wave spacetimes has been analyzed in detail by a number of 
authors [20, 21, 22, 23, 24, 25, 26, 27]. One essential conclusion of this work is that 
nearly all interesting properties of plane wave spacetimes may be deduced from the 
properties of 2 x 2 matrices E = E(^i) satisfying the differential equation 

E = HE. (9) 

This is a "generalized oscillator equation" with — H acting like a matrix of squared 
frequencies. Eq. (9) arises when solving for geodesies or Jacobi fields in plane wave 
spacetimes. Bitensors such as Synge's function and the parallel propagator may be 
written explicitly in terms of its solutions. The same is also true for a plane wave's 
Killing vectors. 

It is convenient to write all possible matrices E in terms of two particular 
solutions. Fix any§ G M and define A{-^Uo) and B(-,iio) to be solutions to (9) 

§ The notation Uo is used here because this wih later be interpreted as a li coordinate associated with 
some observer. Similarly, Us is often interpreted below as a w coordinate associated with a source 
(not necessarily when that source is observable at any observation point with u = Uo)- Later, the 
notation Ue is used to denote a u coordinate associated with the emission of light from a source. 
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(where derivatives are applied to the first arguments of A and B) with the initial 
conditions 

lim A{us^Uo) = lim 9(i)B(ii5, i^o) = 5, (10) 

Ug ^Uq Ug ^Uq 

lim 9(i)A(ii5,'Uo) = lim 'B{us,Uo) = 0. (11) 

Ug yuo Ug yuo 

Here, the notation 9(i)A indicates a partial derivative with respect to the first 
argument of A. We assume for simplicity that H is a matrix of piece wise-continuous 
functions and that A and B are at least (and piecewise-C^) in both of their 
arguments. Example expressions for A and B are discussed in Sects. 5 and 6. 

Given any two solutions Ei and E2 to (9), it is easily verified that their Wronskian 
is conserved: 

E]'E2 - E]'E2 = constant. (12) 

Here, denotes a matrix transpose. Applying this formula with Ei — ^ A and E2 ^ B 
shows that 

ATa(i)B-a(i)ATB = ^. (13) 

Neither A nor B are necessarily symmetric matrices. Nevertheless, (12) and (13) may 
be used to show that the products 

ATa(i)A, a(i)AA-^ BT5(i)B, 5(i)BB-\ BA^, B-^A (14) 

are symmetric wherever they exist [23]. Letting Ei B(-,iio) and E2 B(-,iXs) in 
(12) shows that 

B{us,Uo) = -B^{uo, Us) (15) 

for all Us^Uo G M. This is essentially Etherington's reciprocity law [1, 3, 38, 39]. A 
similar calculation may be used to show that 

A(iXs, Uo) = -5(1) AT(ix^, Us) (16) 

as well. 

It is sometimes useful to consider partial derivatives 9(2) with respect to the second 
arguments of A and B. The resulting matrices remain solutions to (9). Comparing 
initial conditions shows that 

S(2) A(^l5, Uo) = -B{Us, Uo)H{Uo), S(2)B(^l5, Uo) = -A{us,Uo). (17) 

A may therefore be derived from B. The opposite is also true wherever detH 7^ 0. 

Geometrically, B corresponds to the transverse coordinate components of a 
Jacobi propagator describing the evolution of deviation vectors along geodesies passing 
between different pairs of points [23]. B is also related to image distortion. Up to 
an overall time dilation factor, it translates small differences in image position on an 
observer's sky to spatial deviations from a fiducial source point. In the language of 
[11], B is proportional to the Jacobi map. The symmetric matrix 9(i)BB~^ plays a 
similar role, but translates source separations to emission (rather than observation) 
angles. It is proportional to an object typically referred to as the optical deformation 
matrix. These statements are explained more fully in Sect. 4.3. 
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2.2. Conjugate pairs 

It is often useful when working with plane wave spacetimes to consider hypersurfaces 
of "constant phase." Recalling the interpretation of the u coordinate as a phase, let 

Su^ :={peM: u{p) = Uo} (18) 

denote such a hypersurface. Given two points Ps G Su^ and Po G Su^ (with Ug ^ Uq), 
define the multiplicity or "index" of these points to be [20] 

I(jps,Po) := 2 - rankB(^is,^io). (19) 

We also set I{po^Ps) = whenever Uq = Us. It follows from (15) that I{ps^Po) = 
I{po^Ps)- For most pairs of points, / = 0. Such pairs are said to be "disconjugate." 
Pairs Ps-,Po satisfying I{ps^Po) 7^ are instead said to be conjugate with multiplicity 
HPs^Po)' Similarly, we call the pairs Su^^Su^ "conjugate hyperplanes" and the pairs 
of real numbers Us^Uq "conjugate phases" when rankB('Us, i/o) < 2. Despite the 
appearance of (19), the index map / : M x M ^ {0, 1, 2} describes phenomena which 
do not depend on any choice of coordinate system. 

Conjugate pairs as described here are closely related to the conjugate points 
commonly considered in differential geometry and optics. In general, distinct points Ps 
and Po on a given geodesic are said to be conjugate if and only if there exist nontrivial 
deviation vectors along that geodesic which vanish at both po and Ps- In plane wave 
spacetimes, this condition reduces to I{ps^Po) > 0. Defining the multiplicity of a 
pair of conjugate points to be the number of linearly independent deviation vectors 
which vanish at those points, that multiplicity is equal to I{ps^Po)- The concept 
of conjugacy associated with / does not, however, require the specification of any 
particular geodesic. It is uniquely defined even for pairs of points connected by multiple 
geodesies or by none. Indeed, these are the only cases where 1^0. 

All strong lensing effects associated with plane wave spacetimes are related to the 
existence of conjugate hyperplanes. It follows from (19) that every pair of conjugate 
phases Ug , Uq satisfies 

det 'B{us^Uo) = det B('Uo, Us) = 0. (20) 

Finding conjugate pairs and their multiplicities may be viewed as a matter of direct 
computation once H is specified. Alternatively, various Sturm-type comparison 
theorems can be used to make general statements regarding the existence and 
separations of conjugate pairs for various classes of plane wave. See, e.g., Chapt. 
XI of [40] for results relating to mathematical problems of this type and [20] for an 
application to "tame" plane wave spacetimes. More specific examples are discussed in 
Sects. 5 and 6 below. 

The qualitative structure of geodesies in plane wave spacetimes is closely related 
to the index /. First note that every disconjugate pair of points Ps-,Po is connected by 
exactly one geodesic. If two points are conjugate, the number of connecting geodesies 
is either zero or infinity. Sizes of geodesically connected regions may be summarized 

by 

dim [(all points geodesically connected to Po) H Su^] -\- I{ps^Po) = ^- (21) 

A similar relation exists for null cones when Ug. 

dim[(all points connected to Po via null geodesies) D + I{ps^Po) = 2. (22) 

The latter result has been referred to as an "index theorem" in [20]. Eq. (21) 
implies that geodesies emanating from Po and intersecting a hyperplane Su^ that 
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Figure 1: Given a preferred point Po oi" a preferred u = constant hyperplane Su^^ 
plane wave spacetimes naturally divide into a number of open four-dimensional regions 
Mn{uo)- These regions are separated from each other by the hyperplanes S^^(^Uo) 
conjugate to Su^- 

is disconjugate to Su^ form a 3-dimensional region. Indeed, these geodesies fill the 
entire hyperplane. More interestingly, geodesies intersecting a conjugate hyperplane 
with multiplicity 1 fill only a 2-dimensional region on that 3-dimensional surface. 
Geodesies intersecting a hyperplane with multiplicity 2 form a line. Similarly, the null 
cone of a point reduces to a 1-dimensional curve on every hyperplane with multiplicity 
1. A null cone intersecting a hyperplane with multiplicity 2 is focused to a single 
point on that hyperplane. These two cases correspond to astigmatic and anastigmatic 
focusing, respectively. Anastigmatic focusing tends to be unstable in the sense that 
perturbations tend to split a single multiplicity 2 phase into two closely-spaced phases 
each with multiplicity 1. 

In many applications, there exists a preferred point Po, or perhaps a preferred 
hyperplane Su^- Po may, for example, represent the position of an observer at a 
particular time. Fixing this point, the set of hyperplanes conjugate to Su^ divides a 
plane wave spacetime into a (possibly infinite) number of open regions JVn{uo). Let 
JVo{uo) denote the largest connected region containing Su^ and excluding any portion 
of a hyperplane conjugate to Su^- If there exists a smallest ri{uo) > Uq conjugate to 
lio, the surface 5'ri(no) clearly contained in the boundary of A/o(iio)- J^i{uo) may 
then be defined as the largest connected region which includes ^'^^(no) a boundary 
and contains points Ps satisfying Ug > ti{uo) and I{ps^Po) = 0. This continues the 
spacetime "above" S^-^(^u^y Similar constructions may be used to define Tniuo) and 
Mn{uo) for values of n other than 1. See Fig. 1. Fixing a particular nonzero integer n 
and real number Uq^ it is not necessary that Tn{uo) exist at all. In general, the domain 
of Tn is an open subset of R. This domain can be empty for some n. 

The geodesic uniqueness result described above can now be reduced to the 
statement that a point Po is connected to another point Ps 7^ Po by exactly one 
geodesic if and only if there exists some n such that Ps G Mniuo)- Two-point tensors 
like Synge's world function, the parallel propagator, and the van Vleck determinant 
may be defined unambiguously throughout "M x (UnArn)." This excludes from MxM 
only a set of measure zero. 
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2.3. Geodesies 

Beyond the quahtative geodesic structure of plane wave spacetimes discussed above, 
it is not difficult to obtain explicit coordinate expressions for all geodesies. Let F C M 
denote some geodesic and 7 : R ^ M an affine parametrization of it. The vector field 
= {d/dvY generating the characteristics of the gravitational wave is Killing, so 
must be constant on P. If this constant vanishes, T is confined to a hypersurface 
of constant phase. Such a geodesic has the form of a (Euclidean) straight line in the 
coordinates ('y,x). 

Geodesies satisfying 7^^^ 7^ are more interesting. In these cases, the affine 
parameter can always be rescaled such that 7^-^^ = — 1- It is then possible to 
identify that parameter with the phase coordinate u. Doing, so u{^{us)) = Us for 
all Us G M. The spatial components 7 := x(7) of any geodesic are fixed everywhere 
once 7o := 7(1^0) and 70 := 7(1^0) have been specified at some fiducial phase Uq- In 
terms of the matrices A and B defined in Sect. 2.1, 

7(1x5) = A{us,Uo)-io + ^{Us,Uo)io' (23) 

The associated v coordinate of T may be efficiently derived using that fact that the 
vector field (2vdy + x'^diY ^ homothety [23]. This implies that 

v{-f{us)) = v{j{uo)) + ^[i^sius -Uo) +7(^s) -7(^5) -7o -70], (24) 

where 

:= -7a7" (25) 

is a constant. 

Ks is closely related to the conserved quantity on T associated with the Killing 
field ^^ . The (unit) 4- velocity tangent to T is related to 7^^ via 

U: = ^. (26) 



It follows that 

-LU^ = ^- (27) 



If the spacetime is nearly flat (so H 0), Ks reduces to a particle's specific energy 
minus its specific momentum in the direction of the gravitational wave. The limit 
—iaU^ (or his 00) may therefore be interpreted as ultrarelativistic motion 
in the direction of the gravitational wave. By contrast, the limit —iaU^ 00 (or 
hvs 0) corresponds to ultrarelativistic motion against the background wave. 

Other Killing fields present in essentially all plane wave spacetimes may be written 

as 

+ (28) 

where = is any 2-vector with the form 

H = A(-, Uo)S{uo) + B(-, Uo)±iuo). (29) 

For each choice of Uq^ there exists a four-parameter family of such vector fields. Each 
of these is associated with a conservation law. Two such conserved quantities may be 
summarized by 

V{uo) := ^= [A^{us,Uo)j{us) - d^i^A^ {us, Uo)j{us)] . (30) 
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Fixing any this 2- vector is conserved in the sense that it is independent of Ug. 
Another conserved 2-vector may be defined by 

C{uo) := ^= \B'^{us,Uo)^{us) - d^i)B'^ {us,Uo)-i{us)\ . (31) 



In the weak-field limit, V corresponds to the specific momentum transverse to the 
gravitational wave. In this same context, C may be interpreted as the conserved 
quantity associated with boosts transverse to the gravitational wave. It constrains 
transverse displacements. Note that both V and C depend on a choice of Uq- This is 
analogous to the choice of origin necessary to define angular momentum in elementary 
mechanics. 

In stationary spacetimes, it is common to discuss various quantities related to 
gravitational lensing in terms of stationary observers (and often stationary sources). 
While plane wave spacetimes are not stationary, there does exist sufficient symmetry 
to define similarly preferred sources and observers. Two geodesies T and V can be 
said to be "instantaneously comoving" u = Uq when 

ns = i^'s. V{uo)=V\uo). (32) 
These conditions imply that the 4- velocities of both geodesies are parallel-transported 
versions of each other on the constant-phase hyperplane Su^- Note, however, that 
geodesies which are comoving at one phase are not necessarily comoving at any other 
phase. 

As implied by (21), bundles of geodesies are strongly focused on conjugate 
hyperplanes. Consider such a hyperplane associated with a phase Tn{uo) conjugate to 
Uo with multiplicity 1. It is clear from (19) that Bn(i^o) := ^{Tn{uo)^Uo) is a matrix 
with rank 1. There therefore exists a unit 2-vector qn(^o) such that 

cC^{Uo)^n{Uo)=0. (33) 

qn(iio) is unique up to sign. Choosing any Pn('^o) orthogonal to qn('^o), the transverse 
spatial coordinates of all geodesies starting at a given point Po focus to the line 

7(rn) = An7o + '^^Pn (34) 
as they pass through S^^(^uq)' Here, w is any real number. All v coordinates may be 
reached on S^^(^uq) by appropriate geodesies. 

Geodesies starting at Po and intersecting a hyperplane S^^(^Uo) with multiplicity 
2 all focus to the single transverse position 

7(rn) = An7o (35) 
as they pass through S^^(^u^y As in the multiplicity 1 case, all values of v may 
be reached by appropriate geodesies. Eqs. (34) and (35) illustrate explicitly how 
conjugate hyperplanes with multiplicities 1 and 2 are associated with astigmatic 
and anastigmatic focusing, respectively. The former case involves focusing in only 
one transverse direction, while the latter case involves simultaneous focusing in both 
directions transverse to the gravitational wave. 

2.4' Distances 

As noted above, all pairs of points not lying on conjugate hyperplanes are connected by 
exactly one geodesic. There is therefore no ambiguity in ascribing geodesic distances 
to these pairs. In particular, Synge's world function 

(j{Ps^Po) •= -(squared geodesic distance between ps and Po) (36) 
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is well-defined whenever its arguments do not lie on conjugate hyperplanes. Plane 
wave spacetimes constitute one of the few examples where a is known essentially in 
closed form: 

Cr{Ps,Po) = ^{Us- Uo) [ - 2{Vs - Vo) + xJ5(i)B(ii5, Uo)B~^{Us, Uo):>Cs 

+ xJB"^(ii5,iio)A(ii5,iio)xo - 2xJB"^(ii5,iio)xs]. (37) 

This is symmetric in its arguments: (j{PsiPo) = (^{PoiPs)- 

The appearance of in (37) indicates that a tends to diverge when its 

arguments approach conjugate hyperplanes. More specifically, suppose that Us ~ 
Tn{uo) and I{Tn{uo),Uo) = 1. Then, 

(^(Ps^Po) - - ^ (—^ — 7^1 [qn(^o) • (xs - An(iio)Xo)]^ (38) 

is an asymptotic approximation for cr if the bracketed term on the right-hand side of 
this equation is nonzero [23]. It follows from (34) that this expression is valid only 
when there does not exist any geodesic passing from po to a point on S^^(^Uo) with the 
same transverse coordinates as pg. The equivalent result if Ps is near a hyperplane 
Sr^(^Uo) with multiplicity 2 is 

a{ps,Po) ^ — 7^ ) |xs - An{uo)yiof. (39) 

2 \Us - Tn{Uo) J 

Again, this is valid only when there does not exist any geodesic passing from Po to a 
point on 5'r^(no) with the same transverse coordinates as Ps. 



3. Image counting 

One of the most basic questions that can be asked regarding a gravitational lens is the 
number of images that it produces of a particular source. Stated somewhat differently, 
how many future-directed null geodesies connect a given timelike curve (the source) 
to a particular spacetime event (the observer at a particular time)? This may be 
answered using the geodesic structure of plane wave spacetimes summarized above. 

To fix the notation, let Po denote a fixed observation event and T the timelike 
worldline of a point source. Assume that V may be parametrized by an everywhere-C^ 
function 7 : R ^ M. The phase coordinate u serves as a useful "quasi-time" [21] for 
plane wave spacetimes ||, so let 7 satisfy u{j{us)) = Us for all Ug G M. The source's 
worldline is not required to be geodesic. 

This section establishes that under generic conditions, an ideal observer at Po 
may see exactly one image of T from each of the "epochs" Mn{uo) described in Sect. 
2.2. Somewhat more precisely, there typically exists exactly one future-directed null 
geodesic from 

:=T^Mn{uo) (40) 

to Po for each n < such that Mn{uo) exists. See Fig. 2. Recall that the boundaries 
of Mn{uo) depend only on the spacetime under consideration, and not at all on the 

II w is not a "true" time function because V aU is null. Generically, there do exist well-behaved 
functions t : M — R where Vat is everywhere timelike [41]. These are time functions in the usual 
sense, although they do not appear to simplify any computations performed here. Incidentally, the 
existence of such time functions implies that plane wave spacetimes are stably causal. Plane waves 
are not, however, globally hyperbolic [42]. Hypersurfaces of constant t are not Cauchy surfaces. 
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Figure 2: Schematic illustration of the gravitational lensing problem. Images 
correspond to future-directed null geodesies from a timelike source V to an observation 
event Po- The arguments of Sect. 3 show that under generic conditions, exactly one 
image is emitted from each region Mniuo) lying before the observation event. The 
dashed lines correspond to constant-ii hyperplanes conjugate to Su^- 



behavior of any particular source. This allows generic bounds to be placed on time 
delays associated with the various images that may be observed. 

Recall from Sect. 2.2 that each point in is connected to Po via exactly one (not 
necessarily null) geodesic. Synge's function (7(7(1^5), Po) is therefore well-defined and 
explicitly given by (37) for all Ug such that 7(1^5) G F^. Consider instead the rescaled 
function 

T^nKUs) := . (41) 

Uo - Us 

If n < 0, the domain of Hn('^s) is equal to all Us such that 7(1^5) G F^. If n = 0, we 
additionally suppose that Ug < Uq for reasons of causality. Images of F^ produced by 
a plane gravitational wave correspond to the zeros of S^. 

Recalling the form (37) for a, it is clear that depends on F as well as the 
matrices A, B, B~^, and 9(i)B. We have assumed in Sect. 2.1 that A and B are at 
least C^. So is 7. The definition of Mn{uo) ensures that det B(iis, iXo) ^ everywhere 
I]^(iis) is defined. B~^ is therefore and is continuous. 

T^n is also monotonic. To see this, note that (9), (13), (14), (37), and (41) may 
be used to show that 

^n = \[- iaT + + (xo - 9(1)6^7) f ] > 0. (42) 

is both continuous and monotonic, so at most one zero can exist for each n. This 
means that at most one image of a source may reach an observer from each epoch 
Afn{uo)- Exactly one such image exists if some u'^u" are known to satisfy 

j:n{u)<0, J:n{u')>0. (43) 

Such bounds are easily established. 

3.1. Lensing between conjugate hyperplanes 

The simplest case to consider is one where Afn{uo) lies "in between" hyperplanes 
conjugate to Su^- Suppose that both Tn{uo) and Tn-i{uo) exist for some n < 0. This 
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is true in Fig. 2 for n = — 1. More generally, (38) and (39) imply that if there are no 
geodesies connecting Po to either ^{jniuo)) or ^{Tn-i{uo))^ 

lim T^ri{us) = — oo, lim T^ri{us) = oo. (44) 

It follows from these limits together with continuity that is surjective on R. 
Since this function is also monotonic, there must exist exactly one emission phase 
Ue G {rn-i{uo) ^ Tniuo)) such that Po and j{ue) G are connected by a null geodesic. 

Under the same assumptions, projectiles moving on timelike geodesies may be 
thrown from ^{us) G F^ to Po only if Ug < u^. Choosing Us—Tn-iiuo) to be sufficiently 
small (but positive), these projectiles can require an arbitrarily large amount of proper 
time to intersect Po- It is somewhat curious that points 7(1^5) G F^ satisfying Ug > Ue 
cannot be connected to Po by any causal geodesic. Such points may, however, be 
reached by suitably accelerated curves which are everywhere causal. 

3.2. The youngest image 

Next, consider the case n = when there exists at least one conjugate hyperplane in 
the observer's past (as occurs in the example illustrated by Fig. 2). The arguments 
given above imply that if ^{t-i{uo)) is geodesically disconnected from po, 

lim So(iis) = —00. (45) 

The other boundary of the domain of Eq occurs at Us = Uq. Here, 

lim_ T^o{us) = - lim_ . (46) 

This limit clearly tends to +00 if 7(1^0) 7^ ^o- Physically, 7(1^0) 7^ implies that Po 
cannot be connected to 7(1^0) by any null geodesic. Assuming that this is true, there 
must exist exactly one Ue G (r_i('Uo), iXo) such that Po and ^{ue) G Fq are connected 
by a null geodesic. 

3.3. The oldest image 

The results of Sects. 3.1 and 3.2 typically suffice to describe the images formed in 
plane wave spacetimes containing an infinite number of conjugate hyperplanes in an 
observer's past. It is, however, important to consider cases where only a finite number 
of conjugate points exist. This occurs, e.g., for finite wavepackets where H has compact 
support. It is also true of linearly polarized vacuum waves that are tame in the sense 
described in [20]. 

If there is at least one conjugate hyperplane in the observer's past, let N denote 
the smallest negative integer such that tn{uo) exists. If there are no conjugate 
hyperplanes in the observer's past, set = 0. The case illustrated in Fig. 2 
corresponds to A" = — 2 if M-2{uo) extends into the infinite past. Regardless, we 
ask whether there exist any future-directed light rays from F^r to Po- Unlike in the 
cases considered above, sources spend an infinite amount of proper time in JVn{uo)- 

Even in flat spacetime, a source that accelerates for an infinitely long time may be 
causally-disconnected from certain observers via a Rindler horizon. Such phenomena 
can be ruled out here by supposing that there exists some finite constant hZrnin > 
such that 

= ^TM\/aV- [|7(^.)|'+7V^.)H(^.)7(^.)] >^min (47) 
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for all Us less than some cutoff. As is clear from (25) and (26), the left-hand side of this 
inequality acts like the square root of a time dilation factor between the coordinate u 
and the source's proper time. Eq. (47) implies that the source's 4-velocity oc 7^ 
satisfies 

< -tau: < -L= (48) 

Y ^min 

sufficiently far in the past. It therefore excludes sources which experience arbitrarily 
large boosts against the background gravitational wave in the distant past. It is 
satisfied by, e.g., sources whose motion is geodesic sufficiently far into the past. 
Assuming that hZrnin exists, it is clear from (42) that may be chosen sufficiently 
small that 

^n{u') < 0. (49) 

Further assuming that j{tn{uo)) is geodesically disconnected from Po (if 7^ 0) or 
that 7(1^0) is not null-separated from Po (if A" = 0), it follows that there exists exactly 
one Ue smaller than tn{uo) (if A" 7^ 0) or Uq (if A" = 0) such that ^{ue) G Fat is 
connected to Po via a future-directed null geodesic. 

3.4' Total image count for generic sources 

The results just described may be summarized as follows: Suppose that there does not 
exist any past-directed geodesic segment from po to F whose endpoints are conjugate 
(in the usual sense). Also assume that the source and observer are not instantaneously 
aligned with the background wave: ^ liuo)- If there exists an "oldest" phase 
conjugate to Uq^ further require that there be some n^^in > such that the source's 
motion is bounded by (48) sufficiently far in the past. 

For each n < 0, these assumptions imply that an observer at Po sees exactly one 
image of F as it appeared in JVn{uo). This provides a strong bound on the possible 
emission times of different images. If J^n{uo) exists for every negative integer A", an 
infinite number of images are formed. If, however, there is some smallest A" < such 
that JVn{uo) exists, lA'l + 1 images appear at Po- Note that these results depend only 
on the waveform H and the phase coordinate Uq associated with the observer. The 
total number of images is the same for all sources satisfying the hypotheses outlined 
above. 

These hypotheses are generic. Recalling (21), the requirement that F exclude any 
points conjugate to Po along a connecting geodesic is equivalent to demanding that 
the source's worldline avoid certain well-behaved one- or two-dimensional subsets of 
the three-dimensional hyperplanes conjugate to Su^- Similarly, the assumption that 
Xo 7^ l{uo) demands only that F avoid a certain line on Su^- The bound (48) on a 
source's asymptotic 4-velocity can fail to hold only for sources which accelerate for 
an infinitely long time. Violating any of these conditions for a fixed observer would 
require that a source's worldline be quite exceptional. Moreover, we now show that 
I A"! + 1 images appear under even broader (but more difficult to state) conditions than 
those just discussed. 

3.5. Non-generic imaging 

Despite the comments made above, our assumptions on the behavior of a source's 
worldline can be violated in certain cases. Suppose, contrary to these assumptions. 
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that there exists at least one point 7e ^ F which is conjugate to Po along a geodesic 
connecting these two points. An infinite number of geodesies then pass between 7e 
and Po. If the connecting geodesies are null, 7e lies on a caustic of the observer's past 
light cone. Continuous images of point sources - "Einstein rings" - are then formed 
at Po (ignoring the associated breakdown of geometric optics). Such cases are not 
considered any further here. If conjugate points between the source and observer are 
associated with non-null geodesies, discrete images of T appear at Po- In these cases, 
the methods used above are easily adapted to find how many images of arrive at Po- 
As already mentioned, there can be no more than one root for each S^. Depending 
on the details of the system, the associated image from Mn{uo) may or may not exist. 

The case considered in Sect. 3.1 where F^ lies in between successive conjugate 
hyperplanes is the simplest to analyze. Assuming that F does not intersect the caustic 
of Po^ no images are formed of F^ if either ^{Tn{uo)) is connected to Po via a timelike 
geodesic or ^{Tn-i{uo)) is connected to Po via a spacelike geodesic. Otherwise, exactly 
one image exists from this region. 

If there exists an oldest conjugate hyperplane S^^(^Uo) described in Sect. 3.3, 
suppose that the source satisfies (47) for some /^min > 0. There are then zero images of 
Ftv if 7{tn{uo)) is connected to Po via a timelike geodesic. There is exactly one image 
if these points are either geodesically disconnected or are connected by a spacelike 
geodesic. 

The last cases to consider concern images of Fq. First suppose that ^ ^{uo). 
If there exists at least one conjugate hyperplane in the observer's past, one image is 
formed of Fq if either ^{t-i{uo)) is geodesically disconnected frompo or it is connected 
by a timelike geodesic. No images are formed if 7(r_i(i^o)) and Po are connected by a 
spacelike geodesic. If ^ l(uo) and there are no conjugate points in the observer's 
past, condition (47) implies that there exists exactly one image of Fq. 

Cases where the source and observer are instantaneously aligned are more 
interesting. Suppose that x^ = 7(1X0)- There then exists one image of Fq with Ue = Uo- 
Recalling that Uo is not in the domain of Sq, it is possible for a second image to be 
emitted from Fq if = somewhere. This may be seen by noting that 

lim_ So = ^(7(^0)) - Vo. (50) 

Us -^Uo 

Two images of Fq can therefore exist when v{'y{uo)) > Vo and x^ = 7(1X0) • 

If a source includes points which are conjugate to the observer (in the ordinary 
sense), there is no simple result for the total number of images formed. Nevertheless, 
it is always possible to say that the total number of images is less than or equal to 
\N\ +2 if a source does not intersect a caustic of the observer's light cone. 

4. Properties of lensed images 

Plane wave spacetimes typically produce multiple images of each source. Even for 
sources whose intrinsic properties remain constant, these images can appear with 
different spectra, brightnesses, etc. We now compute these properties for generic 
configurations satisfying the hypotheses summarized in Sect. 3.4. 

For each image of a timelike worldline F seen at Po, there is an associated null 
geodesic segment connecting Po to an appropriate emission point 7e = "y{ue) G F. In 
terms of Synge's function (36), these points satisfy 

^(7e,Po) =0. (51) 
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First derivatives of a are always tangent to the connecting light ray. In particular, the 
vector 

a._ ^Mle,Po) , . 

at Po points along the geodesic which eventually intersects 7e (and is therefore past- 
directed). The derivative operator here is understood to act on the second argument 
of a. Also note that is normalized such that iar^ = 1. Parallel-transporting to 
the observation point yields 

Uo - Ue 

The derivative operator in this equation is understood to act on the first argument of 
a. Both and may be viewed as (dimensionless) separation vectors between po 
and 7e. 

Eq. (37) and the various identities of Sect. 2.1 may be used to compute the 
explicit coordinate components of and r^. Components transverse to the direction 
of wave propagation are 

Te = B~^{Ue, Uo) [x^ - S(i)BT(iig, iio)7e] , (54) 

where = x(po) and 7e = x(7(iie)) denote the transverse coordinates of the observer 
and source. A similar calculation shows that 

To = B~^{Ue, Uo) [A{Ue, Uo)yio - le] • (55) 

Yo and Tg are related via 

Ye = 5(i)B(iie, Uo)Yo - d(^i)A{Ue, Uo)yio' (56) 

Much of the discussion below considers sources moving on geodesies. In these 
cases, use of (23) shows that 

Yo = B~^{Ue, Uo)A{Ue, Uo)Syio - To- (57) 

Here, 6x0 := x^— 7(1^0) = x^ — 7^. The various identities involving A and B discussed 
in Sect. 2.1 may also be used to reduce the imaging condition (51) to 

Ks{Uo - Ue) = 2{jo • ^^o - ^Vo) - (5xJB"^(iXe, Uo)A{Ue, Uo)5yio- (58) 

Here, 5vo := v{po) —v{^o)' Eq. (58) is a nonlinear relation for the emission "time" Ue 
in terms of the observer's position po and the parameters 7o, 7o, '^(70), i^s describing 
the source's worldline. As discussed in Sect. 3, there can be many solutions to (58). 
These correspond to different images. 

Neither Ue nor depends on the observer's motion. Nevertheless, redshifts 
and angles on the observer's sky do depend on that motion (as is true even in flat 
spacetime). It is often useful to fix this effect by supposing that the observer is 
instantaneously comoving with the source. Following (32), this is taken to mean that 
the unit 4- velocities Ug^U^ of the source and observer on Su^ satisfy 

eaU:{Uo)=eaUS = -^, (59) 



p,K) = u,K) = Uo = ^ = ^. (60) 



Recall from Sect. 2.3 that the "transverse momentum" Vs{uo) is generated by 
contracting Ug{uo) with the two Killing fields equal to X^.^^ at Po and having vanishing 
first derivative at that point. Also note that (60) implicitly defines an instantaneous 
observer velocity = ^/l^U^ normalized (like 7^) such that £aX^ = — 1. 
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4.1. Frequency shifts 

Gravitational lenses typically discussed in astrophysics involve nearly-Newtonian mass 
distributions which may be regarded as approximately stationary (at least on sub- 
cosmological timescales). If both a source and an observer are sufficiently far from 
such a lens, there can be no significant redshift or blueshift from the gravitational field 
of that lens. Roughly speaking, a light ray falling into any stationary gravitational 
potential must climb out of that same potential. This result breaks down if light 
passes through non-stationary regions of spacetime. Indeed, plane wave spacetimes 
may produce images with significant frequency shifts [19]. 

Consider an approximately monochromatic beam of light emitted from 7e and 
received at A future-directed tangent vector oc — to the emitted light ray 
may always be chosen such that 

UJe = -ke • Ue (61) 

is the angular frequency of the light as seen by its source. The frequency uJq of this 
same light ray as measured by an observer at Po is —ko-Uo^ where /c^ is equal to 
/Cg parallel transported from the source to the observer. The observed and emitted 
frequencies are therefore related by 



Uq ' ko Uq ' 




(62) 



Here, := l/(^ • Ugiue))'^ and := l/(^ • ^o)^- The 2- vectors and appearing 
here are determined by the source and observer positions via (37), (51), (54), and 
(55). The resulting frequency shift is valid for all emission points not contained in a 
caustic of Po- 

Now suppose that a source moves on a geodesic and that the observer is 
instantaneously comoving with this geodesic in the sense of (59) and (60). Then 
= l/(^ • Us)'^ doesn't depend on which image is chosen. Eqs. (56), (57), and 
the symmetry of A may be used to rewrite (62) as 

= 1 + (B-T5x,)T ( J (B-T^x,). (63) 

UJe \f^s + |B T(^Xo|^y 

The matrix in parentheses on the right hand side of this equation acts like a 
metric for the "separation" 2-vector 'B~'^ {ue,Uo)[sio — 7('^o)]- If both eigenvalues of 
A{ue,Us)A'^{ue^Us) — S are negative, the source is necessarily redshifted. Conversely, 
sources are always blueshifted when this matrix is positive definite. If AA""" — S has 
both positive and negative eigenvalues, the sign of the frequency difference depends 
on the direction of B~'''(5xo. For special configurations, there is no frequency shift at 
ah. 



4.2. Angles 

Various images formed from a single source appear at different points on an observer's 
sky. Like redshifts, the relative angles between images change depending on an 
observer's 4- velocity U^. The angle between two images arriving at Po with tangents 
and r^^ is 

^^.f,_ i9ab + Uo,aUo,bKr'o' _ 
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Simplifying, 

cosO = 1- ^r-- — ^— — r-. (65) 

This expression is valid for arbitrary source and observer configurations. Specializing 
to geodesic sources and comoving observers, 

2/^, |(B-iA-B'-iA0(5xo|^ 
" ^ " {^s + |B-iA(5x,P)(,,, + |B-iA^(5x,P) ' ^^^^ 

Here, A = A{ue^Uo) and A' = A{u'^^Uo)' It is evident that angles are largely 
controlled by the difference between B~^A at the two emission times. 

Another interesting angle to consider is the observed separation between a 
single image (emitted at 7e) and a generator ^^ of the background gravitational wave. 
For arbitrarily moving source and observer configurations, 

COS?/^ = l pr. (67) 

For observers comoving with geodesic sources. 

This may be used to rewrite the angle between two different images partially in 
terms of the angles ijj and those images make with . Using (66), 

n , , (B-iA(5xo) • (B'-iA^(5xo) . , . 

COS ^ = COS COS H ^ ^ ^ — i,.^, i ^ / r — ^smw^mw . (69) 

|B~-LA()Xo||B'~-LAMxo| 

Similarly, the frequency shift (63) of an individual image may be rewritten as 

CJo ^ i^s csc^(V^/2) 

It is evident from this equation that images which appear highly blueshifted to 
comoving observers must satisfy ip ^ {). 



4.3. Image distortion and magnification 

Thus far, all sources here have been modelled as though they were confined to timelike 
worldlines. Real objects are not pointlike, however. They form extended worldtubes in 
spacetime. Images of such worldtubes form null geodesic congruences which converge 
on po. These images can be significantly distorted by the curvature of spacetime. 
It is simplest to quantify such distortions by first fixing a particular null geodesic Z 
passing between some part of the source and Po- Precisely which geodesic is chosen 
is not important. Z serves only as an origin from which to discuss nearby light rays 
connecting Po to other points in the source. Once this origin has been fixed, the image 
of an extended source may be described entirely using deviation vectors on Z (at least 
for sufficiently small sources). See Fig. 3. 

Deviation vectors (or Jacobi fields) satisfy the geodesic deviation (or Jacobi) 
equation along Z. Letting denote the past-directed null vector tangent to Z and 
obtained by parallel-transporting from po, every deviation vector is a solution 
to 

r'V5(r^Ver) = R\cdr'r%''- (71) 
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Figure 3: Imaging for an extended source. The fiducial light ray Z is indicated together 
with another light ray separated from it by a deviation vector The vector is 
also drawn. This is tangent to Z at the observation point Po- 



This equation is linear, so must depend linearly on initial data. In particular, 
all deviation vectors can be written as linear combinations of appropriate bitensors 
contracted into the initial data^ ^'^('^o) and (^^{uo)- All light rays observed at Po must 
necessarily intersect that point, so it suffices to set £^^{uo) = 0. The first derivative of 
a deviation vector at Po describes an angular deviation between one point of an image 
and the center associated with Z. We therefore consider deviation vectors with the 
form 

e = B''A{;Uo)i^{uo). (72) 
B"- A is known as a Jacobi propagator. It satisfies the Jacobi equation 

r''Vb(r'=VeB%) = R\cdr''r^B\ (73) 
along Z together with the initial conditions 

lim B^a{u,,Uo)=Q, lim /VbB°^(w„ «„) = (5^. (74) 

Us^Uo Us^Uo 

Note that a is a bitensor. It maps vectors at Po into vectors at others points on Z. 
The transverse components of a are 

BaAX'^i)X^j) = (75) 

where B is the matrix defined in Sect. 2.1. Other components of B^a may be deduced 
from the eigenvector relations [23] 

B'^Aius, Uo)r^ = {us - iio)r^, r^5a^(^., ^o) = - ^o)r^, (76) 

5%(^X„ Uo)£^ = {Us - Uo)t\ LB'^AiUs, Uo) = {Us - Uo)iA^ (77) 

All parts of an image must arrive at an observer along null geodesies. Additionally, 
an observer with 4- velocity can only measure angles of vectors orthogonal to U^. 
It therefore suffices to restrict attention to deviation vectors satisfying 

ro ■ iiuo) = Uo ■ = (78) 

^ Capital letters are used in this subsection to denote abstract indices associated with the observation 
point Po- This is done to avoid confusion when writing down two-point tensors such as B"^ a [see 
(72)]. 
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at Po. These constraints restrict all interesting initial data to a two dimensional space. 
The orthonormal vectors 

form a basis for this space at Po- They satisfy 

^(0 ' ^U) = ^^i' ' e(i) = Uo • e(i) = 0. (80) 

A 

(i) 



Parallel-transporting e^.^ to another point on Z yields 



e?,,=-m-''"-2(;-^=^)r". (81) 

which forms a Sachs basis [1] on Z. Initial data appearing in (72) must be of the form 

e^K) = (C/o-ro)a«e(J), (82) 

where a is an unconstrained 2- vector. The factor {Uq ■ ro) > is included here so that 
a is directly related to angles on an observer's sky. A sufficiently small image may be 
described by a suitable set of 2-vectors a representing the angular locations of each 
portion of the image with respect to the center defined by Z. 

Each a may be translated into a physical displacement at the source using (72) 
and (82). First note that for every particular a, (76) implies that r - ^ = throughout 
Z. Indeed, is always a linear combination of the e^.^^ together with r^. Components 
of proportional to are physically irrelevant, so we consider only the Sachs 
components 

:= e(,) • C = [{Uo ■ ro)BaAe1,)ef^)W^\ (83) 
Defining the 2x2 matrix 

^(i)(j)(^s,^o) := {Uo ■ ro)BaA{us,Uo)e'l^.^{us)efj^{uo), (84) 

it is then clear that ^ = Da for any a. D is referred to as the Jacobi matrix or Jacobi 
map [1, 11]. Using (75)-(77), (79), (81), and (84), 

T>{us, Uo) = {Uo • ro)B{us, Uo). (85) 

If a source moves on a geodesic which is instantaneously comoving with the observer, 

Uo-ro = . (86) 



This discussion implies that a portion of an image with angular separation ea 
from the fiducial direction associated with is spatially separated from the fiducial 
emission point 7e G by 

e^{ue) = {Uo ■ ro)B{ue, Uo){ecx). (87) 

The factors of e <C 1 have been introduced here to emphasize that this description 
is valid only for infinitesimal deviations. Regardless, (87) shows that up to the time 
dilation factor {Uo • r^), the matrix B central to all aspects of plane wave geometry 
may be physically interpreted as a transformation converting infinitesimal angles on 
the vertex of a light cone into infinitesimal separations elsewhere on that light cone. 
B depends only on the u coordinates of the source and emission points, and not on 
any other aspects of the physical configuration. It may be computed for all possible 
observer-source pairs directly from the wave profile H. 
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Angles of emission (as opposed to observation) of the various light rays travelling 
from the source to the observer may be found by differentiating (87) and applying the 
appropriate time dilation factor: 

^ = ) %)BK, ^o)a = j a(i)B(..e, Uo)a. 

The last equality here makes use of (62). Applying (87) shows that 

4 = 5(i)BB-i^. (89) 

The symmetric matrix 9(i)BB~^/([/e • ^e) therefore converts spatial locations to 
emission angles within the source (with the constraint that all light rays intersect 
Po). It is referred to as the optical deformation matrix [1, 11]. 

Eq. (87) implies that there is a sense in which circles on the observer's sky 
correspond to ellipses near 7e. This deformation may be parametrized by performing 
a polar decomposition on 'D{ue^Uo): 

Here, R/3 and represent rotation matrices through some angles /3 and x- The ratio 
Dj^/D- is related to the ellipticity of the aforementioned ellipse, x represents the 
angle between the principal axes of that ellipse and the Sachs basis. D± and x ^^'^ 
referred as shape parameters [1, 43]. 

Recalling that polarization vectors are parallel-transported in the geometric optics 
approximation [3], any polarization vector must have Sachs components which are 
constant along Z. In principle, the angle x might therefore be measured by comparing 
the relative "rotation" between an object's observed shape and an appropriate 
polarization angle [1, 44]. For linearly polarized waves where H can be made diagonal 
by an appropriate coordinate choice, X = with respect to this coordinate system 
and the basis (81). It is shown in Sect. 6 that x ^Iso vanishes in a natural way for all 
sufficiently weak wavepackets which are nonzero only for short times. 

Eq. (87) implies that D converts angles at the observer to separations within 
the source. The determinant of D must therefore relate solid angles at Po to physical 
areas near 7e: 

= |detD| = \D+D_\ = {Uo -rof \detB\. (91) 



It follows that 



4ng := V|detD| = {Uo ■ ro)^/\d^\ (92) 

may be interpreted as an "angular diameter distance." Absolute value signs are 
necessary here because detB changes sign after each pass through a conjugate 
hyperplane with multiplicity 1. Physically, such sign changes represent parity 
inversions of the resulting image. Note that dang does not necessarily increase 
monotonically with the age of an image (as computed using the source's proper time). 
Closely related to the angular diameter distance is the luminosity distance 

<^lum := (^o/^e)"^<^ang. (93) 

One factor of oJo/oJe arises here from considering light cones emanating from the source 
instead of the observer. The other factor of uj/oOe is related to the energy change 
associated with frequency shifts. 
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5. Symmetric plane waves 



Now that various optical quantities have been computed for general plane wave 
spacetimes, we consider their application to various special cases. The simplest 
nontrivial plane waves are the symmetric waves. These are locally symmetric in the 
sense that V aRhcJ = 0- It follows from (3) that symmetric plane waves must have 
constant waveforms. Also note that {d/duY is Killing in these examples [as well as 
= [d/dvY^, which is Killing in all plane wave spacetimes]. Particular symmetric 
plane waves may be specified entirely by the (constant) eigenvalues of H. 

Recalling the decomposition (6) of H into h^ and /ly, a coordinate rotation 
may always be used to set /ix =0 for symmetric waves. It is then evident that the 
two eigenvalues of H are given by — h\\. It is always possible to set 

where 

hi := h\\ + h2 := h\\ - h^. (95) 

The weak energy condition implies that h\\ > 0, so at least one eigenvalue of H must 
be negative (implying that at least one of the hi ^2 must be positive). We assume for 
definiteness that /i+ > 0. Then, 

hi>0, hi>\h2\. (96) 

If the vacuum Einstein equation is imposed, h\\ = and /12 = —hi. For conformally- 
fiat geometries representing spacetimes associated with, e.g., pure electromagnetic 
plane waves, /i+ = and h2 = hi. Other cases may be viewed as superpositions of 
gravitational and ("gravito"-) electromagnetic waves. 

All symmetric waves produce an infinite number of images of almost every source. 
It is clear from (94) that these waves are also linearly polarized. The angles x P 
appearing in (90) therefore vanish when considering image deformations with respect 
to the Sachs basis (81). Other lensing properties depend on the sign of /i2- We call 
the case /12 < "gravity-dominated" and the case /12 > "matter-dominated." 



5.1. Gravity -dominated symmetric waves 

Consider symmetric plane wave spacetimes where /12 = h\\ — hj^ < 0. Gravity- 
dominated waves such as these generalize the vacuum waves satisfying hi = — /i2- 
Symmetric vacuum waves arise from, e.g., the Penrose limit of a null geodesic orbiting 
a Schwarzschild black hole on the light ring. 

For any gravity-dominated symmetric wave, the matrices A and B defined in 
Sect. 2.1 are 

A/ \ ( COS h?(Us—Uo) \ /nr7\ 

y cosh |/i2| 2 [Us - Uo) J 

B{us.Uo) = ( hf'sinhUus-Uo) ^ ^ \ (gg) 

\ I/12I 2 smh|/i2h(iXs - iXo) / 

It is clear that detB(-,iXo) has an infinite number of zeros for any choice of Uq. Each 
of these zeros represents a phase conjugate to Uq. There are an infinite number of 
such phases in both the past and future of every observer. The discussion in Sect. 3 
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therefore implies that under generic conditions, an infinite number of images appear 
for almost every source. Explicitly, all conjugate phases are given by 

_ 1 

Tn{uo) = Uo^ nirh^ ^ , (99) 

where n is any nonzero integer. It is evident from (19) that all of these phases have 
multiplicity 1. For any n < and any observation point Po with u{po) = Uq^ exactly 
one image of each source is visible as that source appeared in JVn{uo). This corresponds 

to the region between u = Uo^nnh^ ^ and u = Uo-\-{n—l)7rh-^ ^ . Note that det B(-, 'Uq) 
switches sign on each pass through a conjugate phase. The parity of an image emitted 
from Mn{uo) is therefore opposite to the parity of an image emitted from Mn-i{uo). 

Specializing to cases where the source is a geodesic and the observer is 
instantaneously comoving with that source on Suoi some configurations lead to 
redshifts and others to blueshifts. Using (63), 

u;e ^K, + |B-T5x„|2- ^'^^^ 

An image is therefore redshifted if and only if 

d^im^o = -hi{8xlf + \h2\{5xlf < 0. (101) 
It is blueshifted when (5xJH(5xo > 0. There is no frequency shift at all in cases where 

{5xl/5xlf = \h^/hr\. (102) 

The direction of the frequency shift clearly depends only on the instantaneous 
orientation 5xo/|^Xo| of the source and the observer on Su^- In particular, it does 
not depend on which image is considered. All images of a particular source experience 
the same type of frequency shift. 

Emission times Ue for an observer comoving with a geodesic source may be found 
by solving (58). For gravity-dominated symmetric waves, this equation reduces to 

ns{uo - Ue) = 2(7o • 5yio - Svo) + {Sxl)'^ cot {uo - Ue) 

+ \h2 I ^ {5xlf COth \h2 I * {Uo - Ue). (103) 

1 

If —n ^ 1, it is evident that the image from Mniuo) must satisfy cot hi {uq — Ue) ^ 1. 
Images from the distant past are therefore emitted at phases Ue very nearly conjugate 

to Uo'- 

Ue^rM-^^^. (104) 

Substituting this relation into (68) and (100) shows that very old images cluster near 
£^ on the observer's sky and experience increasingly-negligible frequency shifts: 

V^(x|n|"\ IcJoM - 1| oc |n|"^ (105) 

Old images of slightly extended sources are also highly distorted and demagnified. 
Their angular diameter and luminosity distances both scale like 



<^ang ^ d\um CX |n| 2 exp 



lV\h2/hi\\n\7ry (106) 



Gravity-dominated symmetric waves therefore produce an infinite number of 
exponentially dimming images for almost every source. 
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5.2. Matter- dominated symmetric waves 

Matter-dominated symmetric waves satisfying /12 > act somewhat differently than 
gravity-dominated waves. In these cases, 

A(..,.„)=f / V (107) 

Y COS /ll ('^s ~ '^o) / 

-- - \ 

h^^ ^mhl{us-Uo) ^ \ ^^Qg^ 

h2 sin h2 {ug — Uo) J 

Phases conjugate to Uq occur at Uq + nirh^ ^ and at Uq + ^ 7 where n, are any 

nonzero integers. If ^Jh\|h2 is an irrational number, these two families of phases are 
distinct. Each conjugate pair then has multiplicity 1. If \/h\lh2 is rational, some 
conjugate pairs have multiplicity 2. In the conformally-flat case where h\ = all 
conjugate phases have multiplicity 2. In every other case where yjhxjh^ is rational, 
an infinite number of conjugate phases occur with each multiplicity. Regardless of 
an infinite number of images are formed for almost every source. 

Now consider a luminous source moving on a geodesic. If the source and observer 
are instantaneously comoving in the sense of (59) and (60), frequency shifts associated 
with each image are given by (63). Since 

B-i(AAT - ^)B-T = H (109) 

is negative- definite in this case, all images are redshifted. 

If —n 1, an image originating from Mn{uo) must be emitted just before the 
source intersects S^^(^Uo)- snoh images cluster towards in the observer's sky 
and have negligible frequency shifts. Images emitted near conjugate hyperplanes with 
multiplicity 1 are highly distorted. Images emitted near conjugate hyperplanes with 
multiplicity 2 are not significantly distorted at all. In both cases, however, older 
images are dimmer (although the rate at which this occurs is much slower than for 
gravity-dominated symmetric waves). 

Regardless of the sign of the oldest images formed by symmetric plane wave 
spacetimes depend on the spacetime structure at arbitrarily large transverse distances. 
If the metric is modified so that the wave decays at large distances, only a finite 
number of images discussed here would be unaffected. The oldest images found in 
pure symmetric waves likely do not appear at all in perturbed symmetric waves. 

6. Sandwich waves 

Symmetric plane waves are mathematically simple, but are not reasonable models 
for gravitational radiation emitted from compact sources. More interesting are waves 
where H is nonzero only for a finite time: Sandwich waves. Suppose, in particular, 
that there exists some ii+ > such that H(ii) = for all u ^ [0,ii+]. It follows 
from (3) that such spacetimes are locally fiat whenever u < {) 01 u > Uj^. The curved 
region containing the gravitational wave is effectively sandwiched between the two 
null hyperplanes Sq and . Every timelike curve eventually passes entirely through 
such a wave. 

Before an observer interacts with the wave, spacetime is fiat and optics is trivial. 
The case Uq > U-^ where an observer has already passed through the wave is more 



B{Us,Uo) 
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interesting. In this case, A and B reduce to their flat space forms 

A{Us, Uo) = 8, ^{Us.Uo) = (us - Uo)S (110) 

when Us > U-^. The forms of A and B inside the wave depend on the details of H, 
and will not be discussed here. If Us < 0, however, there always exist four constant 
2x2 matrices ol^ ol^ (3 and /3 such that 

^{Us.Uo) = (a + OLUo) + (/3 + $Uo)us. (Ill) 

Note that the dots on d and $ do not refer to derivatives in this case. They are only 
used as a labelling device. It follows from (17) that A{us^Uo) is independent of Uq. 
Moreover, 

A{us, Uo) = -d - $Us. (112) 

If there were no wave at all, cx = $ = and — d = (3 = 5. 

In general, d and /3 have a simple physical interpretation. If two geodesies are 
comoving and have a transverse separation 5^o when u > Uj^/ii follows from (23) that 
the transverse separation between these geodesies is — d(5xo immediately before they 
interact with the wave at = 0. Similarly, the relative transverse velocity of these 
geodesies is — /35xo when u <{). 

Interpretations for a and (3 are somewhat less direct. Consider two geodesies 
which intersect at some time Uq > U-^ after the wave has passed, but which have 
a relative transverse velocity Jx^ at u = Uq. The difference in transverse velocities 
between these two geodesies is then {/3-\-Uo$)S±o when < 0. Similarly, the difference 
in the transverse positions of these geodesies is (a + Uo(y.)SiCo at = 0. It follows 
that a controls shifts in position that are independent of the time Uo at which the 
two geodesies cross each other. 

In principle, a, d, /3, (3 may all be found by solving (9) if H is known. At first 
glance, this would appear to imply that 4-4 = 16 numbers are required to describe 
observations through a sandwich wave. The actual number of required parameters 
is somewhat less than this. First note that the Wronskian identity (13) implies that 
a.'^ l3 is a symmetric matrix. It also implies that 

a^(3 = S. (113) 

Further simplifications arise by recalling from (14) that BA""" and B'''9(i)B are 
symmetric. It follows that 

adT, f3f3\ (114) 

are symmetric as well. These expressions are completely general, and hold for any 
sandwich wave. They significantly constrain the number of independent parameters 
needed to specify A and B. Equivalently, they limit the number of parameters that 
must be extracted from H. 

It follows from the arguments of Sect. 3 that the number of images of a 
generic source observable in any plane wave spacetime is governed by the number 
of hypersurfaces conjugate to the u = constant hypersurface Su„ containing the 
observation event Po- Continuing to assume that Uq > all phases conjugate to 
Uo must be smaller than U-^. It follows from (20) and (111) that conjugate phases 
occurring before the wave may be found by solving 

det[(a + cxuo) + (/3 + $Uo)r] = (115) 

for all r < 0. This equation is quadratic, so at most two solutions exist. An observer 
ahead of the wave may therefore see at most three images of a source as it appeared 
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behind the wave. There may also be at most one image of a source as it appeared 
ahead of the wave. In principle, any number of images may arise from inside the wave 
[where (111) is not valid] if H is sufficiently large. 

Conjugate phases found by solving (115) clearly depend on the observation time 
Uo- Less obviously, the number of conjugate phases can also depend on Uq- For 
an observer moving on a timelike worldline (where Uq increases monotonically ) , new 
conjugate phases - and therefore new images - sometimes appear at discrete times. 
These images correspond to observation times where (115) momentarily degenerates 
to a linear equation. New images can therefore arise when Uq = Uq and 

det(/3 + /3iio) = 0. (116) 
If det /3 ^ 0, the two solutions to this equation are 



-[Tr/3Tr/3-Tr(/3/3)] ± A/[Tr/3Tr/3 - Tr(/3/3)]2 - 4det/3det/3 

Uo : . (117) 

2det/3 

Only solutions satisfying Uq > Uj^ > {) are physically relevant. When a conjugate 
phase of this type first appears, it satisfies 

lim t{uo) = — oo. (118) 

The associated image therefore provides a picture of the infinitely distant past. 
Furthermore, an infinite amount of the source's history appears to the observer within 
a finite amount of proper time. This implies that new images are highly blueshifted. 

No matter how long an observer waits, no conjugate phase can exceed Uj^. The 
emission time for an associated image might therefore be expected to tend towards a 
constant value as ^ oo. This means that a very large amount of proper time at 
the observer corresponds to only a small amount of proper time at a source. Images 
which appear suddenly and are initially highly blueshifted become highly redshifted 
at late times. 

It is unclear precisely what these types of images imply. To the extent that 
geometric optics remains valid, all observers passing through momentarily see 
almost the entire universe appear infinitely blueshifted as it was in the infinitely distant 
past. Furthermore, (70) implies that all of the universe is briefly compressed into a 
single point on each observer's sky. Of course, such phenomena lie outside the domain 
of geometric optics. They may even lie outside of the realm of test fields propagating 
on a fixed background spacetime. Extreme focusing events like these might indicate 
instabilities inherent in the plane wave geometry itself. It should, however, be noted 
that all of the infinities just alluded to are likely to have finite cutoffs in "realistic" 
plane waves which decay at large transverse distances. 

Weak wavepackets 

One important class of sandwich waves are those that are very weak and last only for 
a short time. In these cases, a, d, /3, and $ may be expanded as integrals involving 
successively higher powers of H. To lowest order in such a scheme, A and B are 
approximately unaffected by the wave while inside of it. Assuming that —Ug^Uo ^ 
the first corrections to this assumption are 

A{us,Uo) ^ S - Us dwli{w), (119) 
Jo 

'B{us,Uo) ^ {ug — Uo)S -\- UsUo / dw'H.{w). (120) 

^0 
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This approximation is consistent with (13) [and therefore (113) as weh]. In terms of 
the matrices appearing in (111), a ~ 0, — d ~ /3 ~ 5, and 



/ d^H(^). (121) 

JO 

Note that Eqs. (119) and (120) should be apphed with care if H involves many 
oscillations of an approximately periodic function. In these cases, the integral of H 
can be very nearly zero. Terms nonlinear in H might then be significant. 

Assuming that (119) and (120) are indeed adequate approximations for A and 
B, a coordinate rotation may always be used to diagonalize /3. There then exist two 
constants 1-Li and I-L2 such that 

(122) 

In this sense, all sufficiently short gravitational plane waves act as though they are 
linearly polarized [so /3 = x = iii (90)]. In terms of the individual wavefunctions 
appearing in (6), 

Hi = / dw[h\\{w) ^ h+{w)], (123) 
Jo 

/ dw[h\\{w) - h+{w)]. (124) 



/" 



The transverse coordinates have also been chosen such that 

dwhx(w) = {). (125) 

/o 

If a wave satisfies the vacuum Einstein equation, /ly = and Hi = — H2. More 
generally, it follows from the weak energy condition that 

>0. (126) 

Now assume that the integral of /i+ is non- negative, which entails only a minimal loss 
of generality. Then, 

Hi>0, ni>\H2V (127) 

We say that a wave is gravity-dominated if < and matter-dominated if > 0. 
These definitions are closely analogous to those used to classify symmetric plane waves 
in Sect. 5. There, a wave was said to be gravity- or matter-dominated depending on 
the sign of the constant h2 = h\\ — hj^ appearing in (94). 

For weak gravity-dominated wavepackets, there can be at most one phase 
conjugate to an observer satisfying Uq > Uj^. If this exists, it evident from (120) 
that 

r-iK) = - ^"° . (128) 
tilUo — 1 

is conjugate to Uq with multiplicity 1. This equation is valid only if t-i{uo) < 0. A 
conjugate phase therefore exists only for observers satisfying 

UiUo > 1. (129) 

Note that HiUo = 1 is the unique physically-relevant solution to (116) in the gravity- 
dominated case. 

Waves that are matter-dominated (so I-L2 > 0) also admit the conjugate phase 
(128) when Uq satisfies (129). In the conformally-fiat case where Hi = H2, this is the 
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Figure 4: Example emission versus observation times in the presence of weak 
wavepackets. The left figure assumes a vacuum wave, while the right figure is a matter- 
dominated wave with I-L2 = 4Hi/5. Dotted lines represent r-i{uo) and r_2('^o)- The 
wave's location is indicated schematically by a thin grey rectangle. The source and 
observer are placed on geodesies assumed to be comoving after the wave has passed. 
In both cases, hig = 1, 1-Li5x\ = l-Lidx^ — 1/2, and l-i\(^^o • (^x^ — ^Vq) = 3. 



only conjugate phase. Unlike in the gravity-dominated case, the multiplicity of r_i 
is equal to 2 for conformally-fiat waves. In all other matter-dominated cases, r_i has 
multiplicity 1 and a second conjugate phase is admitted (also with multiplicity 1) for 
all observers satisfying 

n2Uo > 1. (130) 

This occurs at 

Note that (130) is a more stringent condition than (129). As implied by the notation, 
r-2{uo) < r-i{uo). 

Consider a point source moving on a timelike worldline F in a weak sandwich wave 
spacetime. Such a source appears differently when observed at different times. It is 




Figure 5: Frequency shifts for a vacuum wavepacket. The parameters and dashing 
used here are the same as those in the left panel of Fig. 4. The younger image 
(solid) experiences a temporary blueshift and then a redshift before settling down to 
ujo = oJe- The slight initial blueshift as Uq 0+ is due to the source and observer 
moving towards each other before interacting with the wave. The older image (dashed) 
initially appears with infinite blueshift. It then suffers an ever-increasing redshift. 
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Figure 6: Bending of emission curves for a vacuum wavepacket. All parameters are the 
same as in the left panel of Fig. 4 except that three different choices are made for the 
value of 5vo' Curves for both images are shifted to the right as 5vo is decreased. The 
younger (solid) curves are almost unaffected by the wave if 5vo is sufficiently large. 
For smaller values of 5vo^ both curves are strongly bent by the constraint that they 
can't pass through the dotted curve representing r_i. 

clear that when < 0, exactly one image of F is viewable under generic conditions. 
As time passes, the wave eventually passes through the observer. A second image then 
appears when Uq = > U-^. This image is always emitted before the first. If the 
wave involves a sufficient amount of Ricci curvature^ (from e.g., electromagnetic plane 
waves) and is not conformally-flat, a third image appears when Uq = > H^^. This 
is emitted before the first two images. All images persist indefinitely once they appear. 
Sufficiently far in the future, one image is observed of the source as it appeared after 
interacting with the gravitational wave. All other images predate this interaction. See 
Fig. 4. 

When the second image first appears at Uq = H^^, the new conjugate phase 
r-i{uo) is divergent. At all later times, it is finite. The same is also true for the 
emission times associated with the second image. Almost the entire past history of 
the source is therefore observable within a finite proper time. This implies an infinite 
blueshift. At late times, r_i — H^^. The observed evolution of the source via 
the second image effectively freezes as asymptotes to —H^^ (which predates the 
source's interaction with the wave). Images such as these are highly redshifted, as 
indicated in Fig. 5. Note that a similar transition from infinite blueshift to infinite 
redshift also applies to the third image if it exists. 

Another qualitative feature of the emission times plotted in Fig. 4 is that there is 
a sense in which pairs of images can "switch roles." Consider, e.g., the left panel of that 
figure. At late times, the solid curve (corresponding to the younger image) is perfectly 
linear. Indeed, it remains very nearly linear until HiUq ~ 5. A rapid transition then 
occurs where the dashed curve effectively takes over this linear behavior while the 
solid curve strongly deviates from it. In a sense, the two images reverse their roles. 
This phenomenon occurs one more time (somewhat less sharply) around HiUo ~ 1 
when the second image first appears. It arises essentially because Ue > r_i for the 
younger image and Ue < t_i for the older image. These constraints can cause emission 
curves to bend sharply - with large accompanying frequency shifts - in order to avoid 

+ Distinguishing between different cases based on the Ricci tensor inside a wave requires that the 
approximations leading to (119) and (120) be vahd. It is possible for, e.g., sufficiently strong vacuum 
waves to admit two conjugate phases in the region u < 0. 
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Figure 7: Angular diameter and luminosity distances for a vacuum wavepacket. All 
parameters are the same as those used in the left panel of Fig. 4. 

intersecting r_i. Whether or not this occurs depends on whether the "average" linear 
increase of ever comes near r_i. If it does, this role switching occurs. If not, the 
younger image is barely affected by the gravitational wave at all. This is illustrated 
in Fig. 6, where emission curves for several sources are plotted simultaneously. 

To summarize, images which appear at discrete times briefly appear as bright, 
highly blueshifted "flashes." Indeed, Fig. 7 shows that their luminosity distances go 
to zero. Simultaneously, the angular diameter distance of each new image tends to 
infinity. It is implied by (70) that all highly blueshifted images make must a very small 
angle ip with on the observer's sky. The second (and third) images therefore appear 
aligned with the direction of propagation of the gravitational wave when they first 
appear. This direction could be quite different from the location of the other image(s). 
As time progresses, all images migrate across the observer's sky as illustrated in Fig. 
8. Different images may remain separated from each other by large angles at all times. 

7. Discussion 

Despite their simplicity, plane gravitational waves behave in qualitatively different 
ways from lenses associated with quasi-Newtonian mass distributions. As expected 
from their dynamic nature, plane waves generically shift the observed frequencies of 




2 4 6 8 UiUo 2 4 6 8 ^^u, 



Figure 8: Angles (in degrees) for images produced by a vacuum wavepacket with the 
same parameters as those used in the left panel of Fig. 4. measures the angle 
between both images, while measures the angle between each image and on an 
observer's sky. 
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various images. They may also admit images which appear to move, deform, change 
brightness, and shift color as time progresses. 

More subtle differences relate to the number of images that are produced of a 
given source. For example, even numbers of images can appear generically (which 
has led plane wave spacetimes to be cited [1] as well-behaved examples where the 
odd number theorem [5, 7] does not apply). Some plane waves can even produce an 
infinite number of discrete images. Perhaps most striking of all are the bright flashes 
shown to be produced by generic sandwich waves in Sect. 6. These correspond to 
individual images which appear at discrete times. More typical gravitational lenses 
can produce new images if a source crosses an observer's caustic. Individual images 
then split into two (or vice versa). The flashes produced by sandwich waves are quite 
different. Their appearance does not require that a source pass through an observer's 
caustic. Such images appear individually from the infinitely distant past. Initially, 
they are infinitely bright and infinitely blueshifted points of light appearing in the 
direction of propagation associated with the gravitational wave. 

Many of these effects depend at least partially on the idealization that a plane 
wave extends undiminished to infinitely-large transverse distances. Plane wave 
spacetimes are not asymptotically fiat. Despite being topologically trivial and locally 
well-behaved, they are not even globally hyperbolic: Null geodesies passing between 
appropriately-chosen pairs of points can extend to arbitrarily large transverse distances 
in between those points. It is this property which permits the infinite number of images 
described in Sect. 5 to be produced by symmetric waves. The fiashes described in 
Sect. 6 also depend on the spacetime structure at arbitrarily large distances. This 
structure likely affects the formation of even numbers of images as well. Indeed, the 
usual proofs of the odd number theorem require global hyperbolicity, among other 
assumptions [1, 5] (see, however, [7] for a more general formulation). 

If a spacetime has the geometry of a plane wave only out to some finite transverse 
distance, all results derived here remain valid if the associated images involve light rays 
which never extend sufficiently far to interact with any large-distance modifications. 
The infinite sequence of images formed by a symmetric plane wave would then be 
expected to become finite for spacetimes which are only approximately plane waves. 
Calculations involving the oldest images could no longer be trusted in these cases. 
Similarly, the bright fiashes associated with ideal sandwich waves are likely to be 
somewhat less extreme for waves which decay at infinity. Large brightnesses and large 
blueshifts can still exist, but these will be cut off at some finite maximum. Such 
maxima may, however, remain quite large. 

It is reasonably clear that modifications of the geometry at large distances can 
remove some images. Less obviously, these modification can also introduce new images. 
Consider, for example, the pp-wave spacetimes obtained by substituting 

Hij{u)x'x^ H{u,x.) (132) 

in the metric (1). These generalize the plane wave spacetimes. It has been shown that 
if H grows subquadratically as |x| ^0, the resulting geometries are globally hyperbolic 
[45]. Moreover, every pair of points is connected by at least one geodesic in these cases 
(unlike in pure plane wave spacetimes where the growth of H is precisely quadratic). 
This implies that modifications of the geometry at large transverse distances can 
introduce new null geodesies even between points at small transverse distances. It 
would be interesting to explore these effects in more depth to understand precisely 
how modifications of this sort (or more general ones) alter the lensing properties 
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described here for ideal plane waves. It would also be interesting to better understand 
what the transient flashes of Sect. 6 imply for waves propagating on plane (or almost- 
plane) wave spacetimes. This can likely be facilitated by the Green functions derived 
in [23]. 
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